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1.0  INTRODUCTION 

SDIO  has  considerable  interest  in  developing  efficient  short  wavelength  Ground  Based 
Lasers  (GBLs)  for  defense  against  ballistic  missiles.  To  date,  SDIO  has  pursued  primarily 
two  GBL  laser  technologies  -  excimers  and  FELs.  In  this  proposal,  Science  Research 
Laboratory  (SRL)  identifies  an  alternate  GBL  concept  based  on  electron  beam  controlled 
discharge  pumped  rare  gas  lasers.  These  lasers  have  overall  system  efficiencies  that  are 
potentially  a  factor  of  2  greater  than  excimer  lasers  and  can  potentially  be  scaled  to  the 
multi-  megawatt  average  power  levels  required  for  strategic  defense  applications.  Other 
advantages  of  these  laser  systems  are  listed  below. 

•  These  lasers  have  known  transitions  in  atmospheric  windows.  Therefore  Raman- 
shifting  the  laser  wavelength  to  wavelengths  that  can  propagate  more  easily  through 
the  atmosphere  -  a  process  which  results  in  decreased  laser  system  efficiency  -  is 
not  required.  For  example  the  EMRLD  XeF  laser  at  351  nm  is  2  -  3%  efficient. 
However  incorporating  a  Raman  shift  to  the  required  411  nm  wavelength  results 
in  a  decrease  in  efficiency  by  about  50%. 

•  With  electron  beam  controlled  discharge  pumping  of  rare  gas  lasers,  most  of  the 
power  and  energy  for  exciting  the  laser  transition  comes  from  the  discharge.  This 
can  be  compared  to  excimer  lasers  in  which  all  of  the  power  and  energy  for  the  laser 
comes  from  an  electron  beam.  The  electron  beam  has  to  be  generated  in  a  vacuum 
and  transmitted  through  a  thin  metallic  foil  to  the  high  pressure  laser  medium.  As 
a  result,  a  significant  efficiency  loss  is  experienced  in  transporting  the  electron  beam 
through  the  foil  and  its  support  structure.  The  fact  that  most  of  the  power  for 
pumping  rare  gas  lasers  comes  from  the  discharge  which  is  controlled  by  a  relatively 
low  current  density  electron  beam  alleviates  the  electron  beam  transmission  loss 
as  well  as  foil  cooling  and  foil  life  issues. 

•  No  fuel  burnup  -  Unlike  excimer  lasers,  these  lasers  do  not  depend  on  halogens 
or  other  fuels  that  are  depleted  by  laser  operation.  Hence  these  laser  systems  are 
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capable  of  long  pulse  operation  and  high  duty  factors  without  gas  processing  and 
replenishment.  Long  pulse  lengths  result  in  a  significantly  lower  peak  power  (at 
a  given  average  power)  which  eliminates  atmospheric  Raman  scattering  losses  on 
both  the  uplink  and  downlink.  Further  the  gases  axe  inert  and  nontoxic  and  so 
chemical  scrubbing  and  special  handling  axe  not  necessary. 

•  Laser  action  occurs  between  allowed  atomic  transitions  and  so  the  gain  is  high  and 
optical  absorption  of  the  laser  medium  at  the  laser  wavelength  is  less  important. 
Therefore  these  lasers  can  be  potentially  scaled  to  ultra-high  average  power  (>  100 
MW)  from  a  single  aperture.  This  feature  will  significantly  reduce  the  cost  of  the 
laser  system. 

•  Efficient  laser  transitions  in  rare  gas  mixtures  occur  between  0.6-4  pm.  At  these 
wavelengths  the  turbulence  in  the  atmosphere  has  a  larger  correlation  length  than 
in  the  ultraviolet  thereby  reducing  the  number  of  channels  required  in  the  adaptive 
optics  system.  Further  at  these  longer  wavelengths,  optical  damage  thresholds 
are  higher  and  atmospheric  transmission  efficiency  is  higher.  These  factors,  taken 
together,  more  than  compensate  for  the  increased  diffractive  spreading  of  the  uplink 
laser  beam  in  terms  of  overall  system  cost. 

•  These  lasers  also  have  the  capability  of  delivering  large  specific  energies  (>  50 
J /liter)  which  will  result  in  relatively  compact  and  inexpensive  ground  based  laser 
devices. 

The  key  issues  which  must  be  addressed  before  the  potential  of  these  lasers  can  be 
realized  is  the  stability  of  the  discharge  as  a  function  of  discharge  pump  power  density 
and  pulse  length.  The  best  results  to  date  have  been  obtained  by  Soviet  scientists  at 
the  Lebedev  Institute/1)  They  extracted  80  Joules  in  a  30  ps  pulse  from  a  5  liter  active 
volume.  More  recently  several  groups  in  the  US  have  obtained  efficient  lasing  action  in 
Ar/Xe  mixtures/2,3’4)  The  most  efficient  and  longest  pulse  data  has  been  obtained  by 
the  scientists  at  Sandia  National  Laboratory/2)  By  nuclear  reactor  pumping  they  have 

2 


SCIENCE  RESEARCH  LABORATORY 


measured  power  efficiencies  of  6%  and  millisecond  laser  pulse  lengths.  The  Soviet  electron 
beam  controlled  discharge  laser  result  is  significant  in  so  far  as  the  specific  energy  is  twice 
as  large  as  excimer  lasers  and  the  pulse  length  is  more  than  a  factor  of  ten  longer  than 
that  achievable  with  XeF  and  XeCl  lasers.  Therefore,  the  same  average  power  can  be 
obtained  from  a  laser  which  has  less  than  half  the  active  volume  of  an  excimer  laser.  The 
increased  pulse  length  will  lower  the  peak  intensity  by  two  orders  of  magnitude  thereby 
eliminating  problems  associated  with  stimulated  Raman  scattering  in  the  atmosphere. 
Science  Research  Laboratory  (SRL)  believes  that  the  Soviet  results  can  be  extended  to 
longer  pulse  lengths  (100/xsec)  resulting  in  larger  energy  extracted  per  unit  volume  (50 
J /liter)  and  shorter  wavelengths  in  the  range  between  0.6  to  2.0  /xm  to  minimize  diffractive 
losses.  The  key  to  realizing  these  performance  improvements  in  lasers  based  on  rare  gas 
mixtures  is  to  develop  stable  electron  beam  controlled  discharge  pumping  at  high  pressure. 
This  key  technical  issues  was  addressed  theoretically  in  Phase  I  by  developing  a  detailed 
kinetic  and  discharge  models  for  these  lasers. 

The  relevant  energy  levels  for  xenon  rare  gas  lasers  are  shown  in  Fig.  1.  Also  shown 
in  Fig.  1  are  the  most  efficient  laser  transitions  observed  to  date.  Table  1  gives  a  more 
detailed  list  of  transitions  that  have  lased  together  with  the  transmission  efficiency  through 
a  3  km  standard  atmosphere  for  laser  radiation  at  each  transition.  Transmission  losses  for 
this  path  are  approximately  equivalent  to  the  uplink  loss  for  a  GBL  located  at  sea  level. 
The  transmission  efficiency  predictions  of  Table  1  are  based  on  high  resolution  atmospheric 
propagation  modeling  obtained  with  FASCODE.  All  transitions  listed  have  acceptably  high 
transmission  except  for  the  5d(5/2);j  — ►  6p(5/2)3  transition  at  2.48  /xm.  There  exist  many 
other  shorter  wavelength  laser  transitions  in  neon  and  krypton  rare  gas  mixtures  which  may 
be  efficiently  accessed  at  high  pressure  with  electron  beam  controlled  discharge  pumping. 

2.0  ELECTRON  BEAM  CONTROLLED  DISCHARGE  EXCITATION  OF 
RARE  GAS  LASERS 

There  are  many  combinations  of  rare  gases  such  as  He/Ne,  Kr/Xe  and  Ar/Xe  which 
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ii-M: 


TABLE  I 

TRANSMISSION  OF  VARIOUS  Xe  LASER  LINES  THROUGH 
3  km  OF  ATMOSPHERE  AT  200  ft  ALTITUDE 

(FASCODE) 
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can  lead  to  efficient,  scalable  short  wavelength  laser  operation  at  atmospheric  pressure.  For 
the  remainder  of  this  section,  electron  beam  controlled  discharge  pumping  of  the  Ar/Xe 
laser  will  be  discussed  in  some  detail.  The  reason  Ar/Xe  is  chosen  as  an  example  is  that 
this  combination  of  rare  gases  has  recently  achieved  an  intrinsic  efficiency  of  5%  in  laser 
experiments  at  1.73  /xm  (1>2«3)and  there  exists  experimental  data  on  high  pressure  electron 
beam  controlled  discharges  in  the  Soviet  literature^1)  which  indicates  that  stable  discharge 
operation  can  be  sustained  for  >  30  /xsec  pulse  lengths  at  the  required  pump  power  den¬ 
sity.  More  Recently  various  groups  in  the  United  States  %3)  including  SRL  have  obtained 
efficient  lasing  action  in  Ar/Xe.  From  this  data  one  can  obtain  a  better  understanding  of 
this  laser.  In  fact  scientists  at  SRL  have  developed  a  model  that  successfully  explains  their 
results.  These  results  and  model  will  be  presented  later  in  this  proposal. 

The  upper  levels  of  Ar/Xe  lasers  are  within  the  Xe*  (5d)  manifold  (see  Fig.  1)  which 
lies  from  10  to  11  eV  above  the  ground  state  of  xenon.  The  resulting  IR  laser  photon 
energy  is  between  0.3  -  0.7  eV,  so  if  these  levels  are  created  by  direct  electron  impact  the 
intrinsic  efficiency  would  have  been  much  less  than  the  5  percent  measured  by  the  various 
groups  that  are  investigating  these  lasers.  To  explain  their  high  intrinsic  efficiencies,  the 
Lebedev  group  have  postulated  creation  of  a  low  level  population  of  Xe*  metastable  levels 
with  a  low  current,  high  energy  electron  beam.  The  low  voltage,  high  current  discharge 
then  ionizes  these  states  by  multiple  electron  impact  and  the  upper  laser  level  is  populated 
by  molecular  ion-electron  recombination.  This  process  is  shown  schematically  in  Fig.  2. 
Typically  the  power  deposited  by  the  electron  beam,  which  is  eventually  channeled  to 
provide  the  required  population  of  Xe*,  is  1-2%  of  the  power  deposited  by  the  discharge. 

3.0  PHYSICS  OF  DISCHARGE  PUMPING 

In  this  section  a  brief  discussion  of  the  physics  of  electron  beam  controlled  discharge 
pumping  of  Ar/Xe  lasers  will  be  presented.  The  high  energy  beam  electrons  e/,  deposit 
most  of  their  energy  into  the  formation  of  electron-  ion  pairs  via  the  following  reaction 
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Figure  2:  Schematic  of  the  pertinent  energy  levels  of  Xenon  High  Efficiency  in  rare 

gas/Xenon  lasers  is  achieved  by  recycling  Xe*  several  times  before  decay  to 
the  ground. 
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e*  +  Ar  — ►  Ar+  +  e  +  e* 


(1) 


The  atomic  argon  ions  Ar+  rapidly  form  molecular  ions  Ar*  by  a  three  body  collision 
process 


Ar+  +  2Ar  — *  Ar£  +  Ar  (2) 

The  molecular  ions  create  xenon  atomic  ions  by  charge  transfer  collisions 

Ar+  +  Xe  ->  Xe+  +  2 Ar  (3) 

The  Xe  ions  then  form  molecular  ions  XeJ  which  rapidly  recombine  with  electrons  to  form 
xenon  metastables  Xe* 


Xe*  +  e  -*  Xe*(5 d)  +  Xe  (4) 

which  then  lase  on  the  5 d  — »  6p  transition  at  1.73  /im. 

hv  +  Xe*(5 d)  ->  Xe*(6p)  +  2 hv 

Once  formed,  the  lower  lying  metastables  of  xenon  will  be  efficiently  ionized  by  multi- 
step  electron  impact  in  the  following  set  of  reactions 


e  +  Xe*  — >  Xe**  +  e 

(5) 

e  A  Xe*  ->•  Xe+  +  2e 

e  +  Xe**  Xe***  +  e 

(6) 

e  +  Xe**  -+  Xe*  +  2e 

e  +  Xe***  -v  Xe+  +  2e 

(7) 
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Since  the  spacing  between  the  manifold  of  excited  states  above  the  metastable  state  is 
less  than  1.6  eV,  the  reactions  (5),  (6)  and  (7)  are  driven  by  the  low  energy  core  of  the 
electron  distribution.  The  fact  that  most  of  the  electrons  can  excite  Xe*  to  higher  lying 
levels,  coupled  with  the  large  electron  impact  cross  sections  for  pumping  the  excited  states 
up  the  energy  level  ladder,  insures  that  reactions  (5),  (6)  and  (7)  proceed  efficiently  and 
rapidly. 

The  upper  laser  level  is  formed  through  the  ionic  channel  as  shown  in  Fig.  2.  Once 
Xe+  is  formed  via  the  multistep  process,  it  can  recombine  rapidly  by  three  body  processes 
to  form  heteronuclear  ions  and  xenon  molecular  ions  via  the  following  reactions. 


Xe+  +  2Ar  — *  XeAr+  +  Ar 

(8) 

XeAr+  +  Ie-*  Xe}  +  Ar 

(9) 

The  molecular  ions  then  recombine  with  electrons  to  form  Xe*  as  shown  in  reaction  (4). 

The  laser  discharge  characteristics  for  a  multi-megawatt  Ar:Xe  laser  are  given  in  Table 
2.  The  results  from  the  Lebedev  Institute  suggest  that  the  optimum  gas  mixture  is  99% 
Ar  and  1%  Xe  at  a  total  pressure  of  3  atm.  The  input  power  per  unit  volume  required 
for  efficient  lasing  has  been  estimated  from  the  results  of  the  Lebedev  research  to  be  104 
W/cm3.  So  for  a  discharge  electric  field  of  600  V/cm,  a  discharge  current  density  of  16 
A/cm2  is  required  and  the  corresponding  electron  density  is  3  x  1014  cm-3.  The  laser 
discharge  pulse  length  is  assumed  to  be  100  /is  which  is  three  times  longer  than  that 
obtained  by  the  Soviet  scientists. 

A  Boltzmann  electron  transport  code,  which  calculates  the  steady  state  electron  energy 
distribution  function,  was  used  to  determine  quantitatively  the  efficiency  and  electron 
impact  excitation  rate  constant  for  pumping  Xe*  up  the  energy  level  ladder.  The  predicted 
efficiency  and  excitation  rate  constant  are  plotted  in  Fig.  3,  4  and  5  as  a  function  of  the 
Xe*  density  for  a  99:1  mixture  of  Ar/Xe  at  a  total  pressure  of  one  atmosphere.  From 
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TABLE 2 


DISCHARGE  CHARACTERISTICS 
GAS  MIXTURE 

TOTAL  PRESSURE 

ELECTRIC  FIELD 

CURRENT  DENSITY 

DISCHARGE  PULSE  LENGTH 

SPECIFIC  ENERGY 

DISCHARGE  CURRENT 

DISCHARGE  VOLTAGE 

DISCHARGE  IMPEDANCE 


Ar/Xe  99/1 
3  atm 
600  V/cm 
16  A/cm2 
100  p.sec 
103  J/L 
5.7  MA 
120  kV 
2x102Q 


E-BEAM  CHARACTERISTICS 


E-BEAM  VOLTAGE 

2  MV 

CURRENT  DENSITY 

0.20  A/cm? 

PULSE  LENGTH 

100  |xs 

PULSE  REPETITION  RATE 

30  Hz 

DUTY FACTOR 

10' 2 

POWER  DEPOSITED  IN  FOIL 

12  W/cm2 

E-BEAM  IMPEDANCE 

200  Q 

GUIDE  MAGNETIC  FIELD 

8  KG 

LASER  CHARACTERISTICS 
WAVELENGTH 

SINGLE  PULSE  ENERGY 
SPECIFIC  ENERGY 
PULSE  LENGTH 
AVERAGE  POWER 
NONSATURABLE  ABSORPTION 
SATURATION  FLUX 
SMALL  SIGNAL  GAIN 
INTRINSIC  EFFICIENCY 
ACTIVE  VOLUME 
OPTICAL  COATINGS 

-  SINGLE  PULSE  FLUENCE  INCIDENT 

-  AVERAGE  INTENSITY  INCIDENT 


1.73  urn 
3.6  x  10 6  Joules 
50  J/liter 
100  jxs 
100  MW 
7  x  10  '  £m  ‘ 

44  kW/cm  2 
0.06  cm  ‘ 1 
5% 

1.2  x  2  x  30  m3 

150  J/cm2 
4  kW/cn? 
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EFFICIENCY 


Xe*  (cm ' 3) 


Figure  3:  Curves  showing  the  efficiency  and  rate  constant  for  electron  impact 

excitation  of  Xe*  ->  Xe**  in  a  3  atmosphere  At /Xe  mixture  containing 
99%  Ar  and  1%  Xe.  The  applied  electric  field  is  300  V/cm. 
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EFFICIENCY 


100 


Xe*  (cm ' 3) 


Figure  4:  Curves  showing  the  efficiency  and  rate  constant  for  electron  impact 

excitation  of  Xe*  -»  Xe**  in  a  3  atmosphere  Ar/Xe  mixture  containing 
99%  Ar  and  1%  Xe.  The  applied  electric  field  is  500  V/cm. 
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EFFICIENCY 


Figure  5:  Curves  showing  the  efficiency  and  rate  constant  for  electron  impact 

excitation  of  Xe*  ->  Xe**  in  a  3  atmosphere  Ar/Xe  mixture  containing 
99%  Ar  and  1%  Xe.  The  applied  electric  field  is  700  V/cm. 
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Fig.  6  which  shows  the  drift  velocity  as  a  function  of  electric  field  it  is  clear  that  the 
current  density  increases  very  slightly  over  this  range  of  electric  field.  The  current  density 
is  proportional  to  the  product  of  (nevj).  For  the  proposed  electric  field  of  600  V/cm  or 
200  V /cm-atm.  The  efficiency  of  exciting  Xe*  up  the  ladder  can  be  high  (>30%)  provided 
the  metastable  density  is  .  1014  cm-3.  For  these  reasons,  SRL  proposes  operating  at  the 
larger  specific  electric  field  of  200  V/cm  atm  which  gives  an  electric  field  of  600  V/cm  at 
mixture  pressures  of  3  atmosphere. 

The  high  energy  electron  beam  current  density  should  be  large  enough  to  insure  that 
the  initial  Xe*  density  is  approximately  1014  cm-3  and  the  electron  density  is  ~  3  x  1014 
cm-3. 

4.0  CONSIDERATION  IN  SCALING  TO  HIGH  AVERAGE  POWER 

From  the  laser  characteristics  listed  in  Table  2,  one  sees  that  it  is  possible  to  extract 
50  Joules/liter  from  the  Ar/Xe  laser  if  discharge  stability  can  be  maintained  for  a  100 
fisec  pulse  length.  To  provide  3.6  x  106  J/pulse  output  energy,  requires  7.2  x  104  liters 
of  active  laser  volume.  Assuming  a  gain  length  of  30  m  and  an  electron  beam  height  of 
1.2  m,  a  flow  width  of  2m  would  be  required.  For  such  an  electron  beam  height  and  at  a 
pulse  repetition  rate  of  30  Hz,  a  laser  medium  flow  velocity  of  1.1  x  104  cm/ sec  (M0. 3)  will 
be  adequate  to  insure  that  the  acoustic  and  entropy  waves  generated  by  discharge  energy 
deposited  in  the  laser  medium  will  be  totally  cleared  by  flowing  gas  before  the  next  pulse. 
A  flow  width  of  1.2m  determines  that  the  e-beam  energy  must  be  2MeV  to  insure  that  the 
electron  beam  and  the  discharge  deposit  power  uniformly  in  the  active  medium.  Such  a 
laser  will  be  capable  of  delivering  an  average  power  of  100  MW  from  a  single  device  which 
is  30  m  in  length  with  the  possibility  of  scaling  to  even  higher  average  power. 

Optical  gain  within  the  narrowline  rare  gas  laser  transitions  is  large  and  it  has  been 
verified  that  the  gain  coefficient  is  substantially  larger  than  the  nonsaturable  absorption 
coefficient,  a  hr.  In  fact  using  the  measured  value  of  a  nr  for  these  conditions  given  in 
Table  2  one  can  estimate  that  a^R  L  ~0.2  for  a  30  m  gain  length  and  efficient  length  scaling 
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to  30  meters  should  not  be  an  issue.  The  appropriate  architecture  for  scaling  to  very  high 
power  is  the  expanding  wave  architecture.  Such  a  laser  architecture,  already  demonstrated 
by  SRL  for  excimer  lasers  under  an  SDIO/ONR  program,  can  be  used  to  maintain  high 
efficiency  as  the  laser  length  is  increased  to  provide  higher  average  power.  The  expanding 
wave  laser  architecture  will  also  reduce  the  effects  of  amplified  spontaneous  emission  on 
the  performance  of  large  scale  rare  gas  laser  amplifiers.  Therefore,  length  scaling  of  rare 
gas  laser  power  amplifiers  is  likely  to  be  constrained  by  medium  homogeneity  and  its 
impact  on  laser  beam  quality  and  ultimately  by  optical  damage.  At  a  laser  length  of  30 
m,  a  medium  homogeneity,  A  p/p,  of  2  x  10-5  is  required  to  preserve  laser  beam  quality. 
This  homogeneity  has  already  been  demonstrated  in  small  scale  flow  loops  developed  for 
excimer  lasers.  Ultimately  laser  length  scaling  will  be  constrained  by  the  fluence  and 
average  intensity  which  can  be  processed  by  the  output  window  of  the  amplifier.  For  a 
conventional  amplifier  architecture,  the  average  output  intensity  is  4  KW/cm2  and  the 
single  pulse  fluence  is  150  J/cm2  for  a  100  MW  laser  design.  At  near  IR  wavelengths,  these 
levels  are  within  current  state-of-the-  art.  These  levels  can  be  reduced  by  a  factor  of  2  by 
using  an  expanding  wave  laser  amplifier  architecture. 

5.0  STABILITY  OF  ELECTRON  BEAM  PUMPED  RARE  GAS 
DISCHARGES 

Discharge  methods  for  pumping  rare  gas  lasers  fall  into  three  categories:  self  sustained 
discharges,  UV /X-ray  preionized  discharges  and  electron  beam  controlled  discharges.  Of 
these,  the  electron  beam  controlled  discharge  is  the  most  stable  and  is  the  only  discharge 
method  that  can  be  scaled  to  the  high  average  powers  required  for  strategic  applications. 

The  reason  that  electron  beam  controlled  discharges  are  more  stable  than  other  dis¬ 
charge  methods  is  that  essentially  all  of  the  secondary  electron  production  is  due  to  the 
high  energy  electron  beam.  The  discharge  electric  field  can  then  be  adjusted  to  optimize 
the  laser  efficiency.  In  self-sustained  and  preionized  laser  discharges,  the  applied  electric 
field  has  to  be  sufficiently  large  to  allow  production  of  secondary  electrons  by  direct  ion- 
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ization  of  the  laser  mixture  by  discharge  electrons  and  also  to  allow  efficient  excitation  of 
the  upper  laser  level.  Typically  the  electric  field  required  for  ionization  is  substantially 
larger  than  that  required  for  the  excitation  of  the  upper  laser  level;  in  fact  it  is  usually  suf¬ 
ficiently  large  to  result  in  instability  and  arcing  and  therefore  the  electron  beam  controlled 
discharge,  which  decouples  ionization  and  excitation,  is  the  preferred  approach. 

Rare  gas  discharges  have  been  investigated  in  detail  for  several  years  and  it  has  been 
shown^  conclusively  that  the  stability  of  these  discharges  are  dominated  by  multi  step 
ionization  processes  in  which  the  first  electron  produces  a  rare  gas  metastable  and  a  second 
electron  ionizes  the  metastables.  Cross-sections  for  ionization  of  the  rare  gas  metastables 
by  electron  impact  are  large,  of  order  10 A2,  since  these  metastable  states  are  electronically 
similar  to  the  adjacent  alkalis  in  the  periodic  table.  From  the  previous  discussions  on  the 
physics  of  discharge  pumping,  the  required  metastable  density  was  5  X  1013  cm-3.  So  for 
a  total  pressure  of  3  atm,  the  ratio  of  the  metastable  density  to  the  ground  state  density 
is  10-6.  Boltzmann  code  calculations  predict  that  the  electron  impact  ionization  rate  of 
the  metastables  is  two  orders  of  magnitude  larger  than  the  ground  state  ionization  rate  at 
these  fractional  metastable  densities. 

The  discharge  model  that  has  been  analysed  during  Phase  I  is  shown  schematically  in 
Fig.  7.  There  are  two  excited  states  nm  and  n 2  and  the  ionization  of  Xe.  nm  is  the  xenon 
metastable  states,  n 2  represents  the  6  p  levels  of  Xe,  and  quasis  neutrality  of  the  discharge 
demands  that  rij  and  the  electron  density  nc  are  equal.  The  time  rate  of  charge  of  these 
three  species  are  given  by 

(10) 


C?72  Tl 

—7^  =  Sm - —  —  7l2nenm  +  721nen2  +  Rn2  +  ( cns)nena 

dt  Tm 


dri2 

dt 


S2  +  7I2  nenm  —  721  nen2  —  Rn2  —  vn2  +  an\ 
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Figure  7:  Schematic  showing  the  various  excited  states  and  reactions 
used  in  the  discharge  model. 
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-jj-  =  seb  +  unen2  -  an\  (12) 

Where  ne  is  the  electron  density,  Seb  is  the  electron  beam  ionization  rate,  nm  is  the 
metastable  density,  a  is  the  recombination  rate,  ( av )  is  the  electron  impact  metastable 
production  rate  constant  and  rm  is  the  metastable  lifetime, S'™  and  S2  are  the  direct  ex¬ 
citation  of  metastables  by  the  electron  beams.  712  is  the  electron  excitation  rates  from 
the  metastable  level  Xe*  to  Xe**  and  721  is  the  reverse  process.  R  is  the  total  decay  rate 
(including  collisional  quenching  and  radiative  decay)  of  n2  to  nm  •  u  is  the  ionization  rate 
of  n2.  In  writing  Eqs.  (10),  (11)  and  (12)  it  has  been  assumed  that  this  ionization  rate 
is  dominant  because  it  requires  only  3  eV  to  ionize  n2  while  it  requires  4.5  eV  to  ionize 
Xe*(nm)  and  13.5  eV  to  ionize  the  ground  state  of  Xe.  Hence  most  of  the  electrons  can 
ionize  n2.  This  fact  combined  with  the  fact  that  n2  has  the  largest  ionization  cross-section 
of  the  three  states  results  in  the  ionization  of  n2  being  the  dominant  ionization. 

Equations  (10)  and  (11)  can  be  added  to  give 

4:(nm  +  ”2  +  ne)  =  Sm  +  S2  +  Seb  -  —  +  { cru)nena  (13) 

dt  Tm 

Seb  »  Sm,S2  and  {av)nena,  hence  the  steady  state  density  of  nm,nmo  is  given  by 

Tm(Seb  +  {<71 /)^eono)  (13a) 

By  substituting  this  value  of  nmo  into  Eq.  (11)  one  finds  that  the  steady  state  density  of 
n2  is 


n2  0  = 


Stbirmrieo  +  1)  +  T~7l2^eona(<7l/) 


721  nt0  +  R 

By  substituting  Eqs.  (13)  and  (14)  into  (12)  one  gets 


(14) 


1  + 


^neo(rm7i2ne()  + 1) 


(721^0  +  K) 


+ 


«eo{ 


^Tm(<7t/)na7i2ne 

721  neo  +  R 


a}  = 


(15) 
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Typically  the  electron  mixing  rate  712 neo  is  very  rapid  and  as  a  result  Eq.  (15)  simplifies 
to 

Seb  [1  +  vTmnt0 ]  +  (vTm(crv)na  -  a)n2eo  -  an\0  =  0  (16) 

So  the  steady  state  electron  density  is  given  by 


neo  =  ]uSebTm  +  {( SeluTm  )2  +  4 Seb(cc  -  {ov)naTm)\ 1/2]  / 2 (a  -  ( cv)narm )  (17) 

From  Eq.  (17)  it  is  apparent  that  for  a  steady  state  solution  to  the  discharge  Eqs. 
(10),  (11),  and  (12) 

a  >  ( <rv)naTm  (18) 

As  the  electric  field  increases  oc  decreases  and  (crv)  increases  and  so  condition  (18)  is 
satisfied  below  some  critical  electric  field  Ec.  For  the  conditions  given  in  Table  2  of  this 
proposal  and  using  the  rate  constants  predicted  by  the  Boltzmann  code  Ec  =  900  V/cm 
so  one  must  operate  the  discharge  below  900  V/cm.  Using  the  rates  predicted  by  the 
Boltzmann  code  (see  Figs.  3-5)  in  Eqs.  (13)  and  (14)  one  can  estimate  the  steady  state 
values  nmo,ri2a  and  neo  as  a  function  of  Seb.  These  results  are  plotted  in  Figs.  8  and 
9.  From  these  figures  one  can  see  that  for  an  Set  of  1020  cm-3  sec-1  the  steady  state 
metastable  nmo  and  secondary  electron  densities  neo  are  consistent  with  Table  2. 

Having  obtained  the  steady  state  densities  of  the  metastables  and  electrons,  the  sta¬ 
bility  of  these  values  will  be  discussed  neet  by  performing  a  perturbation  expansion  about 
the  steady  state  values,  i.e 

^m(0  =  ft mo  4"  ^m(0 

^2  (<)  =  N 20  +  ^2  (t) 

ne(t)  =  Neo  +  f>e(t) 
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SCI 


This  expansion  is  substituted  into  Eqs.  (10)  and  (12).  Asstuning  that  voltage  across 
the  discharge  does  not  vary  one  obtains  the  following  first  order  linearized  equations 

iu>6m  =  —712  (nmo^e  +  neoSm) - —  +  72l(^2neo  +  f>en2o)  +  R&2  +  (&v)na8t  (19) 

I'm 

—iu>82  =  7l2(nmo^e  +  «eo7m)  ~  72l(^2neo  +  t>en2o)  ~  R$2  +  2  ocnto8t  (20) 

— iuj6e  —  v82  neo  +  vn20Se  —  2  aneo8e  (21) 

The  (-zw)  on  the  RHS  of  Eqs.  (19), (20)  and  (21)  are  the  result  of  a  Laplace  transform  of 
the  linearized  equations. 

The  characteristic  matrix  for  this  system  is  given  by 

_  X  +  du  <Zi2  °13 

Dcu  =  0-21  x  +  a22  023  (22) 

0  032  x  +  a33  _ 

where 

Oil  =  7l2™eo  +  1  /t 

Oi2  =  — (721neo  +  R) 
a13  =  7l2n10  —  721^20  +  (ov)na 

021  =  — 712  neo 

0-22  =  721^eo  +R  +  Vn  CO  (23) 

<223  =  ~7l2n10  +  721^20  “  2  OC  7leo  +  1/71 20 

<2-32  = 

a33  =  2cme<>  -  un20 

The  characteristic  equation  in  this  case  is 

(r+o33)[(x+on)(x+a22)— ai2°2i]— a32[(a:+an)a23— °2iai3]  =  x3+Ax2+Bx+C  =  0  (24) 
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where 


A  =  CL\i  +  022  +  033 

B  =  ona22  +  011O33  +  022O33  —  0j202i  —  023032 

C  =  013022033  +  032021013  —  0320iia23  —  012033021 

The  conditions  for  all  the  roots  of  Eq.  (24)  to  have  a  negative  real  part  (stable  conditions) 
are 

B]C]D  =  AB  -C  >  0  (25) 

The  coefficients  A,  B,C  and  D  have  been  computed  for  the  three  cases  considered  previ¬ 
ously  and  the  results  are  shown  below: 


ne 

1012 

1013 

1014 

1015 

fl 

2-107 

2.5-107 

6.7-107 

4.8-108 

E=300  V/cm 

B 

2-1013 

2.7-1013 

1.7-1013 

9-1015 

C 

4-1017 

4.4-1018 

8.3T020 

4.5-1021 

D 

4.4-1020 

6.9-1020 

1.1-1022 

4.3-1024 

D 

2.1-107 

3-107 

1.1  -108 

9.4-108 

E=700  V/cm 

B 

2.1-1013 

3.2-1013 

3.3-1014 

2.3-1016 

C 

4.1-1017 

4.6-1018 

7.2-1019 

4.6-1022 

D 

4.6-1020 

9.5-1020 

3.7-1022 

2.1-1025 

The  conclusion  is  that  the  discharge  is  absolutely  stable  under  the  conditions  con¬ 
sidered. 

To  verify  these  conclusions  the  Eqs.  (10),  (11)  and  (12)  have  been  numberically  inter- 
grated  and  the  results  axe  shown  in  Figs.  10,  11  and  12  for  electric  fields  of  300  V/cm, 
700  V/cm  and  1500  V/cm.  In  all  three  cases  5ej  was  assumed  to  be  1020  cm-3  sec-1.  As 
one  can  see  from  Figs.  10  and  11  a  stable  steady  state  is  reached  for  the  lower  electric 
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fields  300  and  700  V/cm.  The  steady  state  values  of  neo,  nmo  and  n2 o  are  consistent  with 
the  predictions  shown  in  Figs.  8  and  9.  When  the  electric  field  is  raised  to  1500  V/cm  a 
steady  state  is  not  reached  this  is  because  the  condition  given  by  Eq.  (18)  is  violated  and 
no  steady  state  solution  exists. 

6.0  Modeling  of  the  1.73/zm  Ar/Xe  Laser 


In  order  to  gain  insight  into  the  important  parameters  governing  operation  of  the 
1.73yum  Ar/Xe  laser  system  a  simplified  rate  equation  model  is  being  developed.  This 
model  is  outlined  in  Fig.  13.  The  upper  laser  level  (Level  2)  is  pumped  via  recombination 
of  electron-ion  pairs.  Collisions  between  electrons  and  Xe  atoms  cause  transitions  at  the 
rate  t~^  back  to  the  continuum,  at  rate  t~%  from  the  upper  laser  level  down  to  the  lower 
laser  level  (superelastic  collisions),  and  collisional  excitation  at  the  rate  Te j1.  Lumped  loss 
rates  from  both  levels  to  lower  lying  energy  levels  are  included  as  rates  r2-1  and  rf1.  For 
convenience,  the  collisional  rates  out  of  the  upper  laser  level  are  expressed  as  ratios  to  the 
collisional  excitation  rate  into  the  upper  level  from  the  lower  level. 


Asym  = 


Tel 
Te  2 


Fn  = 


Tel 
Te  3 


The  rate  equation  set  then  becomes 


^l  =  an\-gA  +  T-i  [Nl  _  jV2(Asym  +  Fn)]  -  t^N2 
at  fit/ 

+  r^Asym  •  n2  -  +  rf1) 


dne 

dt 


m 

Eeff 


—  an\  +  rel1Fn  •  N2 


where  the  small  signal  gain  g  is  defined  by 


(26) 

(27) 

(28) 

(29) 

(30) 


g  =  [JV2  -  NAgi/g*)] 


(31) 


where  AT2  and  Nt  are  the  upper  and  lower  laser  level  population  densities,  respectively,  ne 
is  the  electron  density,  hv  is  the  photon  energy  (  0.72  eV  ),  a,  is  the  stimulated  emission 

28 


SCIENCE  RESEARCH  LABORATORY 


s 


Source  Function  (Recombination) 


Figure  13:  Simplified  ArXe  laser  model 
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cross  section  and  gu  and  gi  axe  the  upper  and  lower  level  degeneracies.  Completing  this 
system  is  the  cavity  optical  flux  equation  given  by 


Fcavi?  + 


1 

2L 


(In  R  -  Cioss) 


+  <f>t 


spon 


(32) 


where  the  factor  Fcav  is  required  since  in  the  experiments  performed  to  date  the  the  electron 
beam  was  masked  to  pump  less  than  the  full  length  of  the  optical  cavity.  The  effects  of 
amplified  spontaneous  emission  were  minimized  by  this  technique.  The  form  of  the  cavity 
flux  equation  used  as  Eq.  (32)  is  valid  if  the  cavity  round  trip  time  (  2L/c  )  is  short 
compared  to  other  significant  times  scales  of  the  system. 

These  equations  may  be  integrated  numerically  using  standard  techniques  (  e.g.,  a 
Gear  method  ).  The  recombination  coefficient  (a)  may  have  a  time  dependance  since  it 
is  a  strong  function  of  the  electron  temperature.  This  may  be  especially  important  in 
the  afterglow  period  of  the  discharge.  The  results  discussed  here  were  obtained  using  a 
constant  value  of  a.  Numeric  difficulties  may  be  avoided  by  normalizing  the  equation  set 
by  appropriate  scale  factors.  Boldface  dimensionless  quantities  axe  defined  by  the  following 
expressions 


where 


and 


N{  —  N i  (df/fc)  71xfcar?  ^ 

(33) 

=  ^•eP'char 

(34) 

<}>  =  $(k-  nckarhu/(Ts ) 

(35) 

t  =  X/  (£  •  Tlchar ) 

(36) 

i-fi1  =  k  •  ne 

(37) 

nckar  =  {P«l«Eefffl2 

(38) 

ce 

II 

“to 
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where  P(t )  is  the  electron  beam  pumping  power  per  unit  volume  (  Po  is  the  peak  pump 
rate  )  and  Eeff  is  the  averaged  energy  expended  per  electron  ion  pair.  The  factor  ncfcor  is 
the  electron  density  reached  in  steady  state  in  the  absence  of  supraelastic  collisions  from 
the  upper  laser  level  (Fn=0).  The  normalized  equation  set  then  becomes 

=  ne2  -  $  [N2  -  N x{gi/gu)]  +  ne  [Nx  -  N2(Asym  +  Fn)] 

ax 

_ -  (39) 

{kTlckarT 2) 

=  $  [N2  -  Ni(gi/g„)]  +  Asym  •  neN2  -  neNi 
ax 

_ _  (40) 

(kTlcharT\) 

~  =  («/*)(/(<)  -  nj  +  Fn  ■  n.N2)  (41) 

5T  =  *  {F“" [Nz  _  +  (2LkncUr)<>nIt-  C'”s)} 

+  TuTT^  ,  (42) 

(*  ncharhu) 

Finally,  an  expression  for  the  spontaneous  emission  rate  into  the  laser  optical  mode  is 
required.  The  form  derived  by  Siegman  (  written  in  terms  of  the  cavity  flux  rather  than 
the  photon  number  )  is 

_  AT2hi/c2cr,  ,  * 

Yspon  —  g  y 

where  V  is  the  volume  occupied  by  the  mode.  Substituting  for  the  fundamental  mode  cross 

sectional  area  of  a  half-symmetric  (plano-concave)  optical  cavity  allows  the  spontaneous 

emission  to  be  written  as  a  function  of  the  mirror  radius  of  curvature  (i?mtr). 

, _ Nihvc^ _ '  (u) 

*spon  2L(LRmiry^\(i-L/Rmiry/*  y 

These  equations  can  then  be  solved  with  an  appropriate  choice  for  the  pump  rate  function 

m- 

The  dimensionless  output  flux  is  given  by 

=  #(1  -  R)/ 2  (45) 
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where  R  is  the  output  coupler  reflectivity.  The  feedback  coupler  is  assumed  to  be  a  high 
reflector  with  the  total  optical  cavity  losses  included  as  Cioaa. 

The  steady  state  solution  of  Eq.  (39-41)  can  provide  information  on  the  small  signal 
gain  and  the  saturation  flux  as  a  function  of  the  mosel  parameters.  The  small  signal  gain 
is  found  by  setting  =  0  and  solving  the  equation  set  algebraically. 

The  saturation  flux  is  defined  as  the  steady  state  flux  which  reduces  the  the  gain  to 
half  of  the  small  signal  (4>  — ►  0)  level.  The  parameters  Asym  and  Fn  must  be  determined 
via  an  iterative  process  using  experimentally  determined  values  for  the  small  signal  gain 
g  and  the  saturation  flux  <f>aat  at  several  pump  powers  as  a  check.  A  Rigrod  analysis 
of  the  time  integrated  pulse  energies  for  three  different  electron  beam  pump  powers  has 
been  performed.  Values  of  the  parameters  rj,  T2,  k,  and  a.  were  estimated  and  fixed  to  the 
following  values. 


Ti  =  40  nsec 

(46) 

r2  =  270  nsec 

(47) 

a  =  2x  10-13m3sec-1 

(48) 

k  =  2  x  10-13m3sec-1 

(49) 

Regions  in  the  Asym-Fn  plane  which  correspond  to  values  of  g  and  <f)aat  within  some 
selected  range  of  the  experimentally  measured  values  may  be  found  by  imposing  one  further 
constraint.  It  was  found  that  using  a  theoretical  estimate  of  the  stimmulated  emission 
cross  section  (~  10-13  cm2)  resulted  in  gains  a  factor  of  102  too  large.  Therefore,  the  small 
signal  gain  at  the  lowest  pump  power  investigated  (2.2  kW/cm3)  was  fixed  at  0.038  cm-1 
by  adjusting  the  cross  section  for  each  choice  of  Asym  and  Fn.  The  results  of  such  a  process 
are  shown  as  Fig.  14  where  regions  which  predict  values  within  10%  of  the  experiment 
are  delineated  for  each  of  the  pump  powers  employed  (2.2,  6.8,  and  16.2  kW/cm3).  It 
has  been  found  that  choosing  Asym=  1.77  (which  would  correspond  to  Te  =  1.2  eV  using 
detailed  balancing  arguments)  and  Fn=  0.8  and  solving  the  time  dependant  equations  Eq. 
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(39-41)  that  the  calculated  pulse  energies  match  the  experiment  well.  Figure  15  shows 
a  comparison  of  calculated  and  measured  pulse  energies  at  a  pump  power  of  2.2  kW. 
While  the  integrated  energies  are  predicted  the  detailed  transient  pulse  behavior  does  not 
yet  match  well.  Further  work  on  the  model  to  improve  its  capability  is  under  way.  An 
acceptable  model  will  permit  rapid  predictions  of  scaling  behavior. 
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COUPLER  REFLECTIVITY  (%) 


Figure  15:  Comparison  of  calculated  and  measured  ArXe  extraction, 
50  cm  gain  length,  10  m/flat  cavity,  and  0.5%  Xe  in  Ar. 
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CONCLUSION 


SRL  has  analysed  both  the  stability  of  the  rare  gas  discharge  and  the  kinetics  of  the 
Ar/Xe  laser.  The  electron  beam  controlled  discharge  is  absolutely  stable  provided  that 
the  electric  field  is  below  a  critical  value.  From  the  analysis  of  existing  laser  data,  it  is 
clear  that  the  electron  kinetics  have  a  significant  effect  on  the  Ar/Xe  laser  .  Because  of 
the  rapid  electron  mixing  rates,  the  rare  gas  lasers  tend  to  operate  less  efficiently  at  higher 
pump  powers. 
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